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ABSTRACT 

The purpose of this work is to analyze and evaluate the 
effect of jamming waveforms on both coherent and noncoherent 
digital communications receivers. Specifically, random processes 
are utilized as jamming models in which it is assumed that the 
jamming waveforms have been produced by a shaping filter driven 
by white Gaussian noise. Such jamming waveforms are then 
assumed to be present at the input of known receiver 
structures (in addition to- the signals and channel noise 
normally present) , and optimum jamming waveform spectra are 
determined for different receiver schemes and modulation 
techniques . 

Graphical results based on numerical analyses are presented 
in order to demonstrate the effect of different jamming 
strategies on receiver performance. In order to quantify 
receiver performance, bit error probabilities are determined for 
binary modulation systems and symbol error probabilities are 
determined for M-ary modulation systems. In each case, the 
error probabilities are functions of signal-to-noise ratio 
(SNR) and j ammer-to-signal ratio (JSR) . Results show that 
it is generally possible to significantly degrade the 
performance of binary as well as M-ary modulation communica- 
tion receivers by introducing suitably chosen jamming waveforms. 
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I . INTRODUCTION 



The theory of statistical signal detection and estimation 
in the presence of additive white Gaussian noise is widely 
described in many textbooks [Refs. 1,2,3]. Signal detectors 
are typically designed and built to either optimize the 
receiver output signal to noise -ratio, or as is the case with 
digital communications receivers, to minimize the error 
probability . 

While it has been demonstrated that receivers designed 
under a white noise interference assumption tend to perform 
reasonably well even when the interference is not white 
[Ref. 4] , the assumption of white noise interference is often 
invalid, especially when the receiver must operate in a jamming 
environment . 

The goal of this thesis is to analyze the vulnerability of 
certain digital communications receivers designed to operate 
in a white noise interference environment, that must operate 
in the presence of jamming also. The mathematical model of 
the jamming utilized is a colored Gaussian noise process 
whose power spectral density is to be shaped in such a manner 
so as to cause a large increase in the receiver probability 
of error. While it is not always possible to solve certain 
spectral shaping optimization problems, it is possible to 
postulate techniques that intuitively achieve efficient 
utilization of the available jammer power. 
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This thesis is divided up as follows. In Chapter II, we 
present results on colored noise interference effects in 
coherent M-ary Phase Shift Keyed (MPSK) receivers, and receiver 
symbol error probability in the presence of noise and jamming 
is derived. In Chapter III we analyze and determine performance 
of a coherent M-ary Frequency Shift Keyed (MFSK) receiver 
operating in the presence of noise and jamming. Chapter IV 
deals with non-coherent Binary Frequency Shift Keyed (BFSK) 
signal detection in the presence of noise and jamming. The 
performance of the well-known quadrature receiver is analyzed 
under dual channel and single channel operation. In Chapter V 
graphical results are presented and discussed, and performance 
comparisons are carried out. The conclusions and interpreta- 
tions of the results obtained are presented in Chapter VI. 
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II. COLORED NOISE INTERFERENCE EFFECTS IN COHERENT 
M-ARY PHASE SHIFT KEYED MODULATION 

A. SIGNAL DETECTION IN THE PRESENCE OF COLORED NOISE 

The system whose performance is to be analyzed is described 
in Fig. 2.1. The structure shown is the optimum receiver for 
recovery of MPSK modulated data", in the presence of additive 
white Gaussian noise. In PSK modulation, the source (or 
modulator) transmits one of M signals s^(t) , where 

i = 1,2,...,M, over a prescribed time interval. Because in 
transmissions and reception these signals are interfered with 
by noise, at the receiver one observes the signal r(t) rather 
than just one of the transmitted signals. Using hypothesis 
testing concepts, we say that under hypotheses H^ , r(t) 
takes on the form 

H.: r(,t) = Ve S.(t) +w(t) +n (t) 

0 _< t £ T , i 

where for M-ary PSK modulation 

S^(t) = /27 t cos ) (2.2) 

1 = 1 , 2 , . . . , M 
k is an integer 

Here W(t) is a sample function of a white Gaussian noise 
process of Power Spectral Density level Nq/ 2 and n^(t) is a 



( 2 . 1 ) 

= 1,2,. ..,M 
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icjurc 2.]. MPFK Hocoiver 




sample function of a colored Gaussian noise process having 
autocorrelation function K^(x) . We assume also w(t) and 

n^(t) are statistically independent random processes. 

The receiver of Fig. 2.1 is, as previously pointed out, an 
optimiim processor (in minimum error probability sense) when 
n^(t) E 0 . The analysis that follows evaluates the effect 
of on the performance of this receiver. Since 

may represent some form of jamming, the error probability 
expression to be derived can be used to determine the vulnera- 
bility of such a receiver to colored noise jamming, or 
conversely, to determine the colored noise spectrum that most 
effectively causes poor or inadequate receiver performance, 
namely, high error probability. 

The signals S^(t) , i = 1,2,3,. ..,M can be shown to 

have cross-correlation coefficients 

T 

p.. ^ / S.(t) S.(t) dt = cos (2.3) 

-LJ 0 ^ J 

i, j = 1,2, ... ,M 

The receiver takes advantage of the fact that we can express 
the S^(t) functions, i = 1,2,...,M, as an exact (rather 
than approximate) expression of a linear combination of two 
functions (t) and (t) . In other words 
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S^(t) 



n=l 



S . ij; ( t) 
in 



i = 1,2,3,. 






(2.4) 



with 



S . 



in 



T 

/ S^(t) iJ;j^(t)dt- n = 1,2; 

i = 1 , 2 , . . . ,M 



(2.5) 



These basis functions >i'j^(t) and 4*2 (t) can be derived via a 
Gramm- Schmitt orthonormalization procedure (or almost by 
inspection in this case). It turns out that <|^j^(t) and 
(which must be orthogonal) are given by 

, V Cos 2iTkt/T 

4^-1 (t) = (^ 

/T/2 



and 



4^2 (t) 



Sin 2Tfkt/T 
/T/2 



where k is an integer. 

It can be easily shown that 



^il 



= Cos 



2tt ( i-1 ) 
M 



i = 1 , 2 , . . . , M 



(2.7) 



( 2 . 8 ) 



and 



^i2 



= - Sin 



2Tr(i-i; 

M 



= 1 , 2 , . . . , M 



(2.9) 
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We define 



= 2Tr(i-l)/M i = 1,2,.. .,M 



( 2 . 10 ) 



and assuming equal prior probabilties , namely, each signal is 
equally likely to be transmitted, the receiver computes 



£! I r^ 1 - 1,2, . . . ,M 

n=l 



( 2 . 11 ) 



and makes decisions based on which value is largest. Thus 
with 



r^ = / r(t) dt n = 1,2 



( 2 . 12 ) 



we have 



T 

=j^Cos ej / r(t) t^ 2 _(t) dt 

T 

+ ^-Sin6^J / r(t)i(j2^^)‘^^ i = l,2,..,M 



(2.13) 



and using simple trigonometric identities. 



= V Cos(6j_ +n) i = 1,2, . 






(2.14 
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Clearly 



V 







1/2 



(2.15) 



where 



T 

V = / ■ r (t) ip. (t) dt , 

^ 0 ■ ^ 



(2.16) 



T 

V = / r(t) i|j-(t) dt 

^0 ^ 



(2.17) 



and 



n = Tan ^ ^ (2.18) 

c 

B. RECEIVER PERFORMANCE 

Since conditioned on any hypothesis i = 1,2,...,M, 

and are Gaussian random variables, we can obtain the 
c s 

statistics of the appropriate random variables, in the follow- 
ing manner. First, we have 



E{V^/H.} 



T _ 

E{ / [/E S . (t) + w(t) 

0 J 



_ 1 

/E / S . (t) ip it) dt 
0 ^ 



+ n^ ( t ) ] ^ ( t ) 

= /E S.j j 



dt} 



= 12 M 

± ^ / • • • / I- 



(2.19) 
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E{V^/H.} 



also 



Var{V^/Hj } 



and 



Var{V^/Hj} 



In Appendix 



T _ 

E{ / [/E S.(t) + w(t) + n^(t) J (t) dt} 
n J o z 



_ T _ 

/E / S.(t) dt = /E S.^ (2.20) 

Q 3 ^ 3^ 

j — 1/2, 3/. ..,M 



T 2 

= E{ [ / [w (t) +n ( t) ] tjy-, ( t) dt] } 

0 Cl 



T T 

= E{ / / [w ( t) +n„ ( t) ] [w ( t) +n^ ( T ) ] i|j-| ( t ) ij;, ( T ) dtdx } 
0 0 ^ c i i 



N T T 

= + / / K (t-x)i|j-| (t)ijj, (x)dt dx (2.21) 

0 0 



T 2 

= E{ [ J [w ( t) +n (t) ] (t) dt] } 

0 Cl 



T T 

= E{ / / [w (t) +n_ ( t) ] [w (x) +n^ (x) ] (t) ij;- (x) dt dx } 

0 0 <_ 1 1 



N T T 

= -^ + j j K (t-x) (t) (x)dt dx (2.22) 

Z 0 0 



A we demonstrate that 
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T T 



/ / K (t-x)ipi (t)iij, {x)dt dx 

0 0 



T T A ? 

= f J K (t-x) ( t) ( x) dt dx = a 

0 0 ^ ^ c 



(A. 7) 



so that and V conditioned on H. have identical variances, 
c s - j 

Observe also that 



E{ [V^-E{V^/Hj }] [Vg-E{Vg/Hj }]/Hj } 



T T 

E{ / [v/(t) +n^ (t) ] t|)-| (t)dt / [ w( t) + n^ ( x) ] ( x) d' 

O' 0 c ^ 



TT N 



= fl 5( t-x) ( t) ( t) dt dx 



00 



T T 



+ ( j K ( t-x) t|;-| (t) ip- (x) dt dx 

0 0 



(2.23) 



We can observe that the first double integral in Eq. 2.23 
is zero, so that 



E{ [V^-E{V^/H. }] [V^-E{Vg/Hj }]/H . } 



tt 



A 2 



= // (t-x)ti^ ( t) 11^2 ( T ) dt dx = <^i 2 



(2.24) 



We 



demonstrate in Appendix A that in general 2 



zero 



so that V and V conditioned on H. may not be uncorrelated, 
c s 3 
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However we are still able to express the joint probability- 

density function of V and V by using the general form [Ref. 

o . s 

5] of an N-dimensional Gaussian random vector X, namely. 



P^(X) 

where 



(2tt) 



N/2|,|l/2 



exp {-i(x 



(2.25) 



= E{X} 



(2.26) 



and 



A 



= E |(X -m^) (X 



(2.27 



In our case, we have a 2-dimensional problem in which (see 
Eqs. 2.19 and 2.20), 



ra 

—X 



/E 



/E S.^j 



j = 1 , 2 , . . . ,M 



2.28: 



and 



• X 



N 



o 2 

r^^c 



'1,2 



1 . 



N 



o 2 

r^"c 



(2.29; 



so that 



-X 



r^o ^ 2, ,2 2 A , 

+^cJ - °1,2 = 



(2.30) 



It is simple to show that 
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(2.31) 



A 



-1 



N „ 

(-f + o2)/a 



-01^2/"^ 



"^1,2'^^ 

N 



Thus 



,V /H.<''o'V«j' 

c s 3 ■’ 







’v - /ES . , 
C jl 




L“ 


V - /ES . o 
L. s D 2 



N 



(f +a2)/A 



-ai_ 2 /a 



-Qi_2/A 
N 

(f ^o^)/Aj 



V - /ES . , 
c jl 



^ V - /ES . T 
•- s 3 2- 



(2.32) 



with j = 1,2,. ..,M . 

Now we need to obtain from this probability density function the 
joint probability density function of V and n conditioned on . 
This type of transformation [Ref. 6] is well known and can be 
used here to obtain 



'■v,n/H. = VP (vcosn,VEinn/n ) 

J C S J 

+ VP,. ,, /„ (-V Cos n,-V Sin n/H . ) , V > 0, 

^c'^s/^j ^ 0 < n < ^. 



(2.33) 



Using the probability density function of V and V (Eqn. 2.32) 

w o 

yields 






V 


r 

1 


VCosn-v'ESj^ 


T 


a -b 


vcosn-v'ESjj^ ' 


1 


0XP i 

^ (20) 


2 


VSinri“/ESj2_ 




“b a 


VSinn-vE£j2 


J 



^ V J 1 

+ -Z^exp 

v''2ttA 



-vcosn-v'ES j . 
-VSinn-/ES . , 



"t 


a -b 






-b a 





-VCosn-vESj^ 

-V£inri“'' ES . _ 
j2 J 



V > 0, (2.33) 

0 < n < ^. 
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and 



N 2 

a = (^ + a^)/A 



(2.34: 



b = 



2 /A 



(2.35) 



This probability density function can be expressed in the form 






^ exp{— j[a (V Cos n + v^ES . ^ 



^ (2tt) 



+ a (V Sin n - /ES . „) 

J ^ 



•2b (V Sin n - /ES . , ) (V Cos n - /ES . ,)] 

J X 



V 



(2tt) 



_ 2 

exp{— (V Cos n + /ES . ^) + a (V Sin n + /ES . 2 ) 



- 2b (V Sin n + /ES . 2 ) ( V Cos n + /ES .-,)]} 

3 ^ ~1 



(2.36: 



which can be simplified somewhat. 

Observe that the exponential of the first term simplifies to 



2 — 2 
a[V +E-2V/E(S.-|Cosn+S.2Sinn)] -2b[V SinnCosn 

J X J z 



- V/E(S., Sinn+S.2Cosn) + ES..S.2: 



(2.37: 



and the exponential of the second term simplifies to 
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2 2 

a[V +E + 2v/E(Sj^ Cos n +Sj 2 Sin n) ] - 2b[V Sin n Cos n 



+ V/E(Sjj^ Sinn + Cosn) + 



(2.38) 



We can now group certain terms together. Observe from Eq. 2.8, 
Eq. 2.9 and Eq. 2.10 that 



Cos n +Sj 2 Sin n = Cos 9 j Cos n - Sin 0j Sin n 



= Cos ( 6 j + n) 



(2.39) 



also 



Sinn + Sj 2 Cos n = Cos 9jSin n “ Sin 6j Cos n 



= Sin (n “ 9 j ) 



(2.40) 



for j = 1,2,...,M, so from Eq. 2.36 we have 






^ exp{-^[a [V^+E-2V/ECos (0 j+n) ] 



'/(2tt) 



2 — 

- 2b [V Sin n Cos n - V/ESin (n- 9j) -E Cos 6 ^ Sin 0 ^] } 



exp{-4[a [V^+E+2V/ECos (0^+n) ] 

/ T~ ^ ^ 

V (2tt) a 



2 — 

- 2b[V SinnCosn+ V/ESin(n~0j)“ECos0jSin0j]} 
for V ^ 0 and 0 ^ n £ • 

S ince 



(2.41) 



Cos(9j +n) = - Cos(9j +n +tt) 



(2.42) 
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and 



Sin(n - ) 

we have 



= -Sin(n-9j+'n’) 



(2.43) 






^ [exp{-^[V^+E-2V/ECos (9 .+n) ] 



\/'(2'rT) ^ A 



2 — 

- 2b [V Sin n Cos n - V/ESin (n-9 ^ ) -ECos 9 ^Sin 9 ^ ] ] } 



1 2 — 2 
+ exp{— j[a[V +E-2V/ECOS (9 j+n+TT) ] — 2b [V Sin n Cos n 



- V/ESin (n-9 j+T7) -E Cos 9 j Sin9j]]}] (2.44) 

for V ^ 0 and 0 £ n £ it . 

It is apparent from the range of v that the two exponential 
terms can be replaced by a single term with v ranging from 0 
to 2tt . 

Thus , we have 



Ptt „ /tt (V,n/H^) = — exp{-4[a [V^+E-2 V/ECos (9^+n) ] 

V/Ti/ti. j j — z ] 

^ v(2tt) A 

2 — 

- 2b[V SinriCos n-V /ESin(n-9j)-ECos9j Sin9j]]} (2.45) 

for V ^ 0 and 0 £ n £ 2tt. The probability density function of 

n conditioned on H. is obtained via integration of P (V,n/H.) , 

J “ , n 3 

namely , 



24 



oo 



(2.46) 



Returning to our decision rule, (Eq. 2.14), recall that we 
decide based on which 

= VCos(0^+ri) i = l,2,...,M (2.14) 

is largest. 

So, if Hj is the true hypothesis, then a correct decision 
is made if 



V Cos (0j+n) > V Cos (0^+n) ; i=l,2,...,M (2.47) 

i ^ j 

Since Cos x is maximum when |xj is minimum, we see that if 
Hj is the true hypothesis, a correct decision is made if 

l0j + nl < l0j_+nl i = 1,2,... ,M (2.48) 

i ^ j 

Now from Eq. 2 . 9 we know that, 

0 j = 2tt ( j -1) /M 

So Eq. 2.54 is satisfied for n in the region 




TT 

M 



n < 




(2.49) 



Thus, the probability of making a correct decision, given that 
Hj is the true hypothesis, Pr{c/Hj}, is given by 
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Pr{c/Hj } 



(2.50) 




-0 



j 



/ 

7T 

M 




(ri/Hj ) dn 



If we make the variable change 



6 = n + 6j 



(2.51) 



Then Eq. 2.50 becomes 



tt/M 

Pr{c/H.} = / P^/„ (B-eVH )d6 

^ -tt/M 3 3 



(2.52) 



Now from Eq. 2.45 and Eq. 2.46 we have 






V 



/ 



(2tt) 



exp{-i[a [V^+E-2V/E Cos (0 .+n) ] 
^ 3 



2b[V Sinn Cos n-V/ESin(n-0j)-E Cos 0j Sin 0j] ] }dV (2.53) 



0 < n < 2 tt 



so that 



P / (3-0.)/H.) = / — ^ exp{-4[a [V^ + E-2V/E Cos 3] 

n/n- J J 0/2 

^ V (2^) A 



- b[V Sin2(3-0j)-2V/ESin(3-20j)-E Sin20j] ] } dV (2.54) 
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and Eq. 2.52 now becomes 



tt/M ” V 12- 

Pr{c/H. } = / / — z=Z= exp{-i(a[V +E-2V/ECOS 6 

^ -tt/M 0 /,o n2. ^ 

V(2tt) a 

2 — 

- b[V Sin 2(B-9j)-2V/ESin(3-26j)-E Sin 29 j] ] } dV d3 (2.55) 

Since the hypotheses have been assumed to be equally likely, 
we have 



Pr {c } 



1 ^ 

M I Pr{c/Hj} 
^ j=l ^ 



(2.56) 



so that 

P = 1 - Pric} 

e 

M 7T /M oo 

= 1 - s I / / — exp{-i[a[V^E 

j=l -./M 0 

- 2V/E Cos 3-b [V^ Sin 2 (3-9^ ) -2V/ESin ( 3-29 ^ ) -E Sin 29^] ] }dV dS 

(2.57) 

Observe that if colored noise is not present, then from Equations 

2 

2.24, 2.30 and 2.35, A = (N^/2) and b = 0, so that Eq . 2.57 
simplifies to the well-known expression for the performance of 
the M-PSK receiver operating in the presence of additive WGN. 

That is. 
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1 



p 



e 




- 2V/E Cos 6+E] ] }dV dS 



(2.58) 



where in Eq. 2.57, the dependence on the index j disappears 
when b = 0. While Eq. 2.57 yields a mathematical result on 
the performance on the M-PSK receiver in the presence of 
WGN and colored noise jamming, its further analysis represents 
a separate project in itself. Not only must Eq. 2.57 be 
optimized for energy constrained jamming but also it must be 
evaluated when the jamming spectrum takes on some simple forms. 
For this reason, no effort has been made to further develop 
the above results . 
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III. COLORED NOISE INTERFERENCE IN COHERENT M-ARY 
FREQUENCY SHIFT KEYED MODULATED SYSTEMS 



A. SIGNAL DETECTION IN THE PRESENCE OF COLORED NOISE 

The structure of the demodulator whose performance is 

to be analyzed is shown in Fig_. 3.1. This receiver is 

known to be optimxim for deciding with minimum probability of 

error, which one of M different signals forming an orthogonal, 

equal energy set received in additive white Gaussian noise was 

actually transmitted. The problem analyzed here, can be 

stated as follows: A waveform r(t), received in the interval 

(0,T), contains one of the M signals, Sj^(t), i = 1,2,...,M, 

with equal probability, as well as white Gaussian noise w(t) 

of Power Spectral Density level N^/2 and colored Gaussian noise 

n (t) having autocorrelation function K (t) . The signals are 
c c 

orthogonal with energy e. That is 



The decision rule used by the receiver, is to choose Sj^(t) 
as the transmitted signal if G^ is a maximum, where 




T 

/ S . (t) S . (t) dt 
0 ^ ^ 




(3.1) 



T 




0 



/ r(t) S^(t) dt 



1 



1,2 



M . 



(3.2) 
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3d 



Figure 3.1. MiFSK receiver 



While this is an optimum test (in minimum probability of error 

sense) in the absence of the colored noise n (t) , the analy- 

c 

sis of the next section is carried out in order to determine 

the effect of n (t) on the receiver performance. Since n (t) 

c c 

will typically be inserted in the channel by an unfriendly 
jammer, it is reasonable to assume that n(t) and n^^(t) are 
statistically independent random processes. 

B. RECEIVER PERFORMANCE 

Since is the output of the ith correlator, and, condi- 
tioned on any hypothesis, is a Gaussian random variable, 
we can obtain the appropriate conditional statistics that allow 
determination of P^, namely the receiver error probability. 



Thus 



T 



E(G./H.} = E{ / [S. (t) + w(t)+n^(t) ]S . (t)dt) 

J n J 



0 



T 



/ S . (t) S . (t)dt 
0 ^ 




(3.3) 



and 



Var{G^/H^} 



T 2 

E{[ / [ w(t) + n (t) ]S . (t)dt] } 

0 ^ J 



tt 



E{// [w(t) + n_(t)] [w(x)+n(T)]S.(t)S.(T) dtdx } 

0 0 ^ J J 



31 



r ) T T N 

Var]G./H.[ = // 6(t-f) + K (t-x) ] S . (t) S • (x) dtdx 

/ D 1) 0 0 2 ^ ^ 

NT T T 

= ^ / S . (t)S . (t)dt + // K (t-x)S . (t)S. (x)dtdx 

^ Q 3 3 00 ^ ^ ^ 

N TT 

= -^ e + II K (-t-x)S . (t)S . (x)dt dx (3.4) 

2 00 ^ ^ 



Define 



„ Tt 

a . = // K (t-x) S . (t) S . (x) dt dx 

c,j oO ^ ^ ^ 



(3.5) 



so that 



Var{Gj/H^} 



N 

o 

2 



£ 



+ 



2 

a 

Of 3 



(3.6) 



Observe furthermore that 



E{ [Gj-E{Gj/H^}] [Gj^-E{Gj^/H^}]/H^} 



T T 

= E{ / [w ( t) 4n _ ( t) ] S . ( t) dt / [w ( X ) +n_ ( X ) ] S, ( X ) dx } 
0 ^ ^ 0 c X 



T T N T T 

= // ( t-x ) S. ( t) S, (x ) dtdx +r / K ( t-x ) S. ( t) S ( x ) dtdx 

OO 2 3 k ^0 3 k 
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E{ [Gj-E{Gj/H^}] [Gj^-E{Gj^/H^}]/H^} 

NT T T 

= ^ / S.{t)S^{t)dt + // K (t-T)S . (t)S, (T)dt dx 

^ 0 ^ 00 ^ 

N TT 

= e 6.^ + If K (t-T)S . (t)S, (x)dt dx (3.7) 

As can be seen from Eq. 3.7, due to the presence of the colored 
noise, the random variables {Gj/H^} are not uncorrelated. 
However we will show that for MFSK with signal frequencies 
that are sufficiently separated, the integral 

TT 

// K (t-x) S . (t) (x) dt dx (3.8) 

00 3 

vanishes for j k, so that the random variables are indeed 
uncorrelated . 

Thus, conditioned on Hj^, the G^ are statistically indepen- 
dent. Assume now that Sj^(t) is transmitted and G^ = x. Then 
the conditional probability of a correct decision, Pr{c/Hj^, 

G^ = x} becomes 
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P{c/H^,G^ =x} 



P{Gj^ <x,...,G^_j^ 




<x/H^ 



G. 

1 



= x} 



M 

n P{G, < x/H. ,G. = x} 
k=l ^ ^ ^ 

k^i 




2’’‘T®k'"°c,k' 



2^2 

exp{-Y /2 (-^e+a^^j^) }dy 



(3.9) 



Introducing a change of variable, 



Z 




(3.10) 



we have 



P{c/H^,G^ =x} 



M 

n 

lc=l 

k?^i 





1 - z ^/2 , 

e ' dz 

/2tt 



(3.11) 



Now, since 



E{Gj_/H^} = £ 



(3.3) 



and 



Var{G./H. } 
1 1 




+ a 



c , k 



(3.6) 



we have that 
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oo 



1 



P{c/Hj_} = 



/ P{c/K^,G^=x} 



27T{- 



N 

o 

T 



? ^ 7 

-exp{-(x-£) /2(r^+a . ) }dx 

^ Cj / X 



c,i 



so that using Eq. 3.11, we obtain 



P{c/H^} = 



oo M 

/ n / 

5 k=l 
k?^i 



xAZ-fcta^^k 



1 -z /2 , 
v/2T ^ 



p IN ^ 

expj-(x-£) /2(-^z+o . ) |dx (3.12) 

Assume now for convenience that M is odd, and express the M-afy 
FSK signals as 



I N 2 



S^(t) = A Cos (OJ^ + (i- (M+1) /2) Aoj) t 0 < t < T (3.13) 

i = 1 , 2 , . . . , M 

so that 



T -jut oo -jut 

/ S.(t)e dt = / SAt)p(t)e 

0 -oo 

= F{S|(t)p(t)} i = 1,2,. ..,M (3.14) 

where 
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p(t) 



1 0 < t < T 

P(w) 

0 otherwise 



-jwT/2 

Te 



Sin (jjT/2 
ooT/2 



(3.15) 



and S|(t) is just Sj^(t) with -<» £ t Thus 



00 

F{S|(t)p(t)} = ^ / S|('co-A)P(A)dA 

— oo 



(3.16) 



where with 



O) . 



1 



(i 



(M+1) 

2 



) Aoo 



(3.17) 



we have 



S!(oj) = ttA [ 6 (oj- (oj +00 • ) ) + 6 (oj+ (cio +00 • ) ) ] (3.18) 

J_ ^ J. CJ ^ 

for i = 1,2,. ..,M . 



Thus 



F{S!(t)p(t)} = j ttA[6 (oo-v- (oo + 00 - ) ) +6 (ou-v+ (oo + 00 - ) ) ]Te 

1 ZTT ' Cl Cl 



-jcoT/2 



Sin vT/2 
vT/2 



dv 



2 



- j (to-oo^-oo^) T/2 Sin (oj-oo^-w^) T/2 
( 00-00^-00^) T/2 



+ e 



— j ( 00 + 00 ^+ 00 ^ ) T/ 2 Sin ( oo+oo^+oo ) T/ 2 
(oo+oo^+oo^) T/2 



(3.19) 



For convenience, let 
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L(.oj) 



-jcoT/2 Sin aiT/2 
® coT/2 



(3.20) 



So that Eq. 3.79 becomes 



F{S!(t)p(t)} = F{S^(t)} = S^(w) = 

AT 

= ~[L(a)-co^~C0j^) + L(co+u)^+(jj^) ] (3.21) 

i = 1 , 2 , . . . , M 



Let us examine now the correlation coefficient p . . , namelv 

ID 




T 

/ S . (t) S . (t)dt 
0 ^ 



T 

/ ACos (oj^+oj . )■ tACos (w + . ) t dt 

' Cl c j 



,2 

A T 

2 



Sin (co^-o)j ) T 

(o) . -0) . ) T 
ID 



+ 



Sin (2cj +to . +C 0 . ) T 
c 1 j 

( 2o)^+a) . +C0 . ) T 
c 1 3 



(3.22) 



If we assume that to^T >> it, then the second term in Eq. 3.22 
vanishes and we have 




T 

/ S. (t)S . (t)dt 

0 J 



2 

A T Sin (i- j ) AcoT 
2 (i-j)AwT 



(3.23) 



In order to have orthogonal signals we need at least AcoT = 
or equivalently Ato = tt/T. Normally, we will have 



Aoj = Rtt/T 



(3.24) 
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where k is a large integer, so that p^j - 0 for i 7^ j . Thus, 
from Eq. 3.23 and Eq. 3.24, 



13 



/ S^(t)S (t)dt = 



£ = A T/2 i = j 
0 i 7 ^ j 

2 



From Eq. 3.5 it appears however that the term a 

c , 1 

independent of i. Nevertheless Eq. 3.12 becomes 



(3.25: 



P{c/H^} = 



00 M 

/ n 

0 k=l 

l^i 




dz 



2 c , 1 



7^7 

exp{-(x-e) /2(-:^a 

2 c , 1 



(3.26: 



Let 



x-e 




Then 



(3.27) 



X 




(3.28) 



So that Eq. 3.26 becomes 
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P{c/H^} = 



00 M 

/ n 

0 k=l 

Mi 



''2 c,i ''2 c,k 



1 ^-z /2 , 

e ^ dz 

/2tt 



1 -n /2 ^ 

X e ^ dn 



(3.29) 



Finally , 



M 



P{c} = 



M I P{c/H.} 
i=l 



(3.30) 



or equivalently 



M 00 M 

^ I / n erfc 



• n , -1 

1=1 -oo k=l 

k?^i 




(3.31) 



We must focus on the a , factor. Observe from Eq. 3.7 that 

c , k 



E{ [Gj-E{Gj/H^}] [G^-E{G^/H^}]/H^1 



N 

^r- £ 6 . 



T T 



_ - + // K (t-T) S . (t) S. (t) dt dx 

2 jk c : k 



(3.7) 



and the second term becomes 
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(3.32) 



T T X oo 

// K (t-T)S.(t)S, (T)dtdT = // K (t-T)S'{t)S;(x)dtdT 

00 J ^ c 3 K 



2tt 



00 

/ 

— oo 



S^(0J)S^ {-cj)Sj^(w) 



do) 



It has been shown in Appendix S that S', (-u)) and S/ (w) are 

3 ^ 

essentially frequency disjoint, therefore Eq. 3.32 is zero 
for j 7 ^ k. For j = k, we have (using Eq. 3.21) 



02 

c,k 



oo 

^ (03) |Sj;.(w) Tdoj 



= ^ (“HL (03-03^-03^) + 1(03+03^+03^^) I ^ d03 ( 3.33) 



If we define 



Ik = ^ / s ( 03 ) I L(03-03^“03j^) + L(o3+03^+03j^) | ^ do3 (3.34) 

—CO 

k = 1,2,. ..,M 

2 

then, with e = A T/2 (Eq. 3.25) 

2 

0 , = el, k=l,2,...,M (3.35) 

C ^ K K 

Thus from Eq. 3.31 
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P(=t = I I 



M 00 M 

/ n erfc 
i=l -00 k=l 

Mi 




- r ?/2 



dn 



(3.36) 



so that 



= 1 - 



M 

- y 

M ‘■ 



oo N 



i=l -oo 



/ n erfc^ 



k=l 

Mi 




(3.37) 



Observe that 






?£+eI . 

2 1 



V 



N 

T-*^\ 



1 



N 

^+el. 

2 X 



V eN 

± 

V 



+£l. 



£N 



/ SNR 



: SNR 



V 



l+2Ij^ SNR 

I+2I! SNR 
1 



3.38) 



Where £/N = 



Then 



I . 

1 



SNR and I ! = 

1 £ 



ith channel JSR. 



M 00 M 

= 1-^1 / I erfc, 

“ i=l -oo k=l 

Mi 



/ SNR 
^ Vo.S+r SNR 



/l+2lj^ SNR 

'I+2I! SNR 
1 



e ^ dn (3. 39) 

/2tt 



Consider now the following colored noise power spectral 
density. 
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(3.40) 



M 



S (oj) 
c 



2ttK y [6((jj+ai +(jj . ) + 6(cu-u -lo.)] 

x=l Cl Cl 



Thus, the colored noise consists of equally weighted "tones" 
at the signal frequencies. Therefore, Eq. 3.34 becomes 



T ” M 2 

I, — I 2t\K y [6(uH-a) +0).) |L((jo— cu "to,,) +L(oj+aj +Wi,) I doj 

K „ -I Cl Cl c K c K ' 

-00 x=l 



M 



2 2 
= TK ^ [ |L(-2a)^-C0j^-a)^) +L(cOj^-oj^) | + |l(co^-coj^) + L ( 2(jj^+c0j^+u^ ) | ] (3.41) 



i=l 



Since is typically large, we can justify the statement that 
the terms involving 2co^ are negligible small, so that. 



M 

= TK[ ^ |L(Uj^-(u^) 1^ + iL(aj^-t^) 1^] 

i=l 



M 



M /Sin(oj.-co, )T/2\ 2 
1 K 



2IK n i = 2TO J I (c..-co,)T/2 

1=1 1=1 \ 1 k ^ 



M /Sin 

= ™ i, (ti- 



M /Sin(i-k)A(i)T/2 \2 



L-k) AcoT/2 



(3.42) 



With AtoT/2 = miT where m is large, we have 



I, = 2TK for i = k 
k 



(3.43: 



We can impose a constraint that 
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M 

^ [6 (orf-o) +to- ) + 6 ((jj-(jj “oj . ) 
i=l = ^ 






= 2KM 



(3.44) 



Then 



and 



Furthermore 



and since 



K = 



P ./2M 
nj 



(3.45: 



P . 
2T ^ 



2M 



TP 



n: 



M 



(3.46) 






_k 

£ 



TP 



nj 



Me 



(3.43: 



TP • = jammer energy and 



n: 

£ 



= signal energy. 



this implies that TP •/£ = JSR. We have therefore that 

n3 

Eq. 3.39 becomes 



P 

e 




M “ p 
I / 1-erf c 
i=l 



•k 




SNR 

5+JSR- SNR/M 



M-1 






dn 



= 1 



oo r f 

-J [l-erfcWn-^y 



M-1 



SNR 

0 . 5+JSR* St4R/M 



- n ^/2 , 
e ' dn 



/2tt 



(3.47) 
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Observe that for JSR = 0, Eq. 3.47 is identical to the well- 
known formula for the performance of the receiver of Fig. 3.1 
under MFSK modulation. 
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IV. NON-COHERENT BINARY FREQUENCY SHIFT KEYED 

SIGNAL DETECTION IN THE PRESENCE OF COLORED NOISE 

A. THE QUADRATURE RECEIVER, EQUIVALENT FORMS AND RECEIVER 

PERFORT4ANCE IN THE PRESENCE OF WHITE GAUSSIAN NOISE 

In this section, a short presentation of the basic princi- 
ples of statistical communication theory that lead to the design 
of the well-known quadrature receiver is undertaken. Basic 
results that are useful in the sequel are presented only, since 
the details have been worked out in numerous textbooks (see 
[Ref. 7] for example) . 

Consider a binary digital communication system model in 
which one of two signals, S^(t) or Sj^(t), with energy E^ and 
E^, respectively, is received in the time interval (0,T). At 
the receiver, white Gaussian noise with zero mean and spectral 
density N^/2 is added to the signal. The actual received 
signal r(t) takes on one of the two forms, namely 

r(t) = /Ej^S^(t) + n(t) 0_<t<_T, i = 0,l (4.1) 

The likelihood ratio test which operates on r(t) in order to 
choose which one of the two hypotheses is believed to be 
the true one, namely 

H^: r(t) = vE“ Sj_(t)+n(t), 0£t<T, (4.2) 

i = 0,1 
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IS 



T 2 

exp{- ^ j [r(t) S^(t)] }dt 

A(r(t)) = — I Y (4.3) 

exp{- / [r(t) -/E. S (t) ] }dt 
0 ° ° 

where y is a threshold whose v^lue depends on the decision 
criteria used. This test can be applied to any communication 
problem involving transmission of known signals S^(t) and Sj^(t). 
One such example is the well-known BFSK modulation scheme. 

One problem of interest, which is a slight modification of 
BFSK modulation problems, involves signals 



/Ej_ S^(t) = A Sin (o)^t+(f.^) i = 0,1 (4.4) 

0 _< t < T 

where the phases (})^, i = 0,1 are statistically independent 
random variables, uniformly distributed over the interval 
(0,2tt) , and the amplitudes A are known and equal. It turns 
out that the test specified by Eq. 4.3 can be modified to 
account for the random phases 4)^ by using conditional proba- 
bility densities. 

The details of the procedure have been worked out in Reference 
8. It can be shown that when the signals are given by Eq. 4.4, 
the test of Eq. 4.3 becomes 



A(r (t) ) 



I^(2Aq^/No> , 



(4.5) 
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where 



2 



' T 12 

/ r(t)Sin 0 ), tdt + 
0 ^ 



T 2 

/ r(t)Cos 0 ), tdt 
0 ^ 



k = 0,1 



(4.6) 



and 1 q ( * ) is the modified Bessel function, defined by 



Io(x) 



00 



I 

n=0 



2n 



X 



2^^(n!)^ 



2it 




xCos ( 0+a) j 
e da 



(4.7) 



For minimum error probability decision criterion, the decision 
rule of Eq. 4.5 assuming equal prior probability of trans- 
mitting S^(t) or , is to choose if 

I (2Aq,/N ) > I (2Aq /N ) (4.8) 

o ^1 o — o ^o o 

or equivalently, to choose if 




Otherwise H is chosen. (Observe that I (x) is a monotonically 
o o 

increasing function.) 

The receiver structure that implements the test of Eq. 

4.8 is shown in Fig. 4.1. Another (equivalent) form of the 
receiver of Fig. 4.1 is shown in Fig. 4.2, involving a combina- 
tion of matched filters and the envelope detectors. The 
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f(t) 




Figure 4.1. Quadrature receiver 




Figure 4.2. 



Matched f ilter equivalent to quadrature 
receiver 
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receiver of Fig. 4.2 is completely equivalent to the receiver 
of Fig. 4.1. 

The evaluation of the performance of the receiver has been 
worked out in Reference 9 and is given by 

, -E/2N 

Pg = 2 ® 

2 

where E = A T/2 is the average signal energy. If we now 
define the signal to noise ratio, (SNR) as 

SNR = E/Nq 
we obtain the simple result 

i exp{- SNR/2} (4.10) 

B. RECEIVER PERFORMANCE IN THE PRESENCE OF COLORED NOISE 

The receiver presented in Section A is optimum in minimum 
probability of error sense when operating in a white Gaussian 
noise interference environment. In this section we analyze the 
vulnerability (probability of error) of the quadrature receiver 
in the presence of an additional additive noise that is modeled 
as colored and Gaussian, having autocorrelation function 
( X ) . (We denote n^(t) as this additional colored noise) . 

The problem can then be restated as follows. Under 
hypotheses H^, i = 0,1, r(t) takes on the form 
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(4.11) 



r(t) = S^(t) + w(t) + n^(t) i = 0,1 

0 < t < T 



where 



S^(t) = /E Sj(t) = A Sin (co^t + (})^) i = 0,1 



(4.12; 



In order to determine the effect of n (t) on the receiver 

c 

probability of error, we evaluate the statistics of the random 
2 

variables k = 0,1, where, as defined by Eq. 4.6, 



2 



T 


2 


T- 


/ r ( t) Sin CO, tdt 


+ 


/ r(t)Cos 00 , tdt 


_ 0 ^ _ 




. 0 ^ J 



(4.6) 



k = 0,1 



Thus, conditioned on H^, i = 0,1 



/ [S^(t) +w(t) +n^(t)]Sin t dt 



+ 



/ [S . (t) +w(t) +n (t)]Sin oj, t dt 
/I c k 



1 2 



A 2 2 

= X. , + Y. , 

1 , k 1, k 



i = 0,1 k = 0,1 



(4.13: 



Observe first that the integral 
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T T 

/ S • (t) Sin 0 ). t dt = / A Sin (oj ■ t + (J) . ) Sin oj, t dt 

0 0 1 1 K 



AT 

2 



"Sin ( CO T/2 
(cOi-cO]^) T/2 



Cos [ (co^“C 0 j^) T/2 






Sin (co j^+coj^) T/2 
(co^+co^) T/2 



Cos [ ( cOj^+-co^) T/2 + (j)^] 



i = 0,1 



k = 0,1 



(4.14) 



If we now assume that 



(co,~co^)T = 2m7T and (co,+to )T = 2Zn 
1 o 1 o 



(4.15) 



we have that Sin(co^ +coj^)T/2 = 0, for i = 0,1, k = 0,1. Thus 



T 

/ S^(t)Sin (Oj^tdt 



2 



Sin(co^-ojj^)T/2 

(c0i-(Ok)T/2 



Cos [ (‘^j_“<\) T/2 + 6^ ^ 



i = 0,1 k = 0,1 



(4.16) 



where 



1 if i = k 

^ ik ” 

^ 0 if i / k 

By arguments similar to the above. 



(4.17) 



T 

/ S^(t)Cos (Oj^tdt 



AT 

2 



Sin(co^-Ci^)T/2 



Sin[ (co^-oj^)T/2 +4>^j 6^ ^ 



i = 0,1 k = 0,1 



(4.18) 
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are Gaussian 



Now conditioned on (J) . , i = 0,1, the X. , and Y. , 

X X ^ K X ^ r 

random variables, so it is possible to obtain the conditional 
probability density function of qj^, k = 0,1. Thus 



E{X. , /h. ({).}= / S. (t)Sin ( 0 , t dt i = 0,1 k = 0,1 (4.19) 

X / K/ X / X Q X K 



and 



r T 



Var{X^^jyH^,(f)j^} = / (w(t)Hn^(t) )Sin (^tdt 



rp r£ 



Jf EK(t)w(T) +n (t) n (t) j-Sin oxtSin oiidt di 
00 c C K K 



T N 2 TT 

/ Sin 0) tdt + // K (t-t)Sin co, t Sin co, xdt dx 

0 2 K 00 ^ K X. 



= 0^+0^, i = 0,1, k = 0,1 
w c,k ' ' 



(4.20) 



where, assuming that >> 1 



o A N T - NT 

2 A o 1-^.2 , ~ o 

a.. = ^ j Sin 0 ), tdt = —j— 



w 



0 



k = 0,1 



:4. 21) 



and 



2 A 
^ c , k 



T T 

= // K (t-x)Sin w, t Sin 00 , xdt dx k = 0,1 

00 K K 



(4.22) 
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Similarly 



E 







T 

/ S. (t) Coso),tdt i = 0,1 k = 0,1 (4.23) 

0 ^ ^ 



and 



Var(YiyH.,*.) 



T N TT 

/ ^Cds^oxtdt + // K (t-x)Cos (jo,t Cos oxxdt dx 

o2 K 00^ ^ ^ 



A 



2^2 
0 + 0 , 
w c,k 



k = 0,1 i = 0,1 



(4.24) 



since it can be demonstrated that 



T T T T 

// K (t-T)Sin ooutsinco, dtdi = // K (t-i)Cos oj, t Cos Co, tdtdi (4.25) 

O'O C K 0(D ^ K K 

k = 0,1 



Finally 






( T T 

e| / [w(t) + r^(t) ]Sin lOj^tdt / [ w(t) - f (t) ] Cos og^t d-; 



T N T T 

o ^ . 



/ -T- Sin 00 , t cos Co, tdt + // K (t-x)Sin ux t Cos oo, xdt dx 

0 2 K K 00 K K 

i = 0,1 k = 0,1 



(4.26) 
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It can be shown that these two integrals are zero so that 
X. , and Y. . are conditionally uncorrelated, and therefore 

1 f K 1 

independent since they are Gaussian random variables. Now 
define 



I 




A 





0,1 



and 



o 



2 

k 




a 



2 

c ,k 




(4.27) 



(4.28) 



so that the conditional density functions for k = 0,1 are, 



(q'+A' )' 
o oo 






oo 



2a 



o 



2o 



o 



o 



g2 

o 



|u(q;) 



(4.29: 



where u(*) is the unit step function, and 



A ' 

oo 



E^{X /H ,4) } + E^{Y /H ,(J) } 
0 , 0 ^ 0^0 0 , 0 ^ 0^0 



Using Eq. 4.16 and Eq. 4.18, we have 



'oo 



/AT ^ , 2 ,AT ^ . ,2 

(-^ Cos i))^) + —Sin 4) ) 

2 o 2 o 



(^)2 

k 2 



(4.30) 



Also 



p(qyH^,0^) = 



exp { — > I 



2o 



o 



2a 



o 



K^lo 



o 



u(q;) 



(4.31) 
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where 



4 ,^} = 0 (4.32) 



due to the result of Eqs. 4.16, 4.17 and 4.18 for i k. 
Therefore 



P(qyH^,4,j_) = 



exp{- qo/ 2 o^}u(q^) 
2 a 

o 



(4.33) 



Furthermore 



P(qi/^o''^o^ = ~T ^ 



2 a^ 






2o^ 






Wv 



(4.34) 



where again due to Eqs. 4.16, 4.17 and 4.18, 



' 0,1 






= 0 (4.35) 



so that 



P(q^/Ho^4>o) 



2 a 



exp{-q^/ 2 a^ }u (qp 



(4.36) 



Finally 



P(q|/H^,4)^) = 



■ exp 



2a 



-(qj_+ A{^) 

2 a? 






11 

2 



u(q[) 



(4.37) 
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(4.38) 



Hi 






(^)2 
^ 2 ’ 



We now have the statistical information needed to compute the 
probability of receiver error P^. Assuming that each hypothe- 
sis has equal prior probability, we have 



= i P{q^ -q >0/H^} + ip{q^ -q^ <0/Hp 



(4.39: 



Observe that 






where 



P{qi >qo/H^,q^ = p} = / P(q^/H^)dq^ 

P 



(4.41) 



Since the conditional probabilities functions are not 
dependent on the individual phases, that is 



P(q^/H^) = P (q^/H^,(^)^) i = 0,1 



(4.42) 



and 



P(qo/H^) = P(qyH^,(J^) i = 0,1 



(4.43) 
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we can rewrite the conditional probability functions in 
the following form 



P(q^/Ho,V 



= 2q P 






= —exp ]- 



((M-/2)^+q^) 



2a 



2 - 



1 ^°( 



qQ(M-/2) 



u(qo) 



P (qo/H]_ = 2q^P (q^H^ , ) 



= -|exp{- q^2a^}u(q^) 






2 2 

= — exp {- q^/2a^}u(qj_) 

^1 






% 

— exp 
^1 



((AT/2)^ +q^) 



q^(AT/2) 



u(q. 



Thus 



oo oo 



p^q^-qQ > 0/Ho) = / [ / 

-oo p ^ 



(4.44) 



(4.45) 



(4.46) 



(4.47) 



(4.48) 
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Similarly 



00 

P{q^-q^<0/H^} = / P|q^ <q /H^,q^ = p^P (p/H )dp 

-00 ^ 



where 



P 

p{qi < qo/^i''5o = = I P(q3_/Hj_) dq^ 

— 00 ' 



so that 



P{q^-q^ < 0/H^} 



oo 

/ 

— oo 



p 

/ P(q^/H^)dq^ 



oo 




(p/H^)dp 



Using now Eq. 4.44 and Eq. 4.46 we have 



pfqi%>o/Ho} = / 



2 2 

/ ^ exp {-q,/2o, }u(q^)dq, 

P a. 



P (P/H )dp 



/ exp{p^/2a^}-^exp 

^o 



2 2 



I /Pl^^/il\u(p)d 



2a 



o I 2 I 

^o f 



For convenience, let e = AT/2 and , recalling u(p) = 1, p > 
Eq . 4.52 becomes 



(4.49) 



(4.50) 



(4.51) 



(4.52) 



0 
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OO ^ X 

P{qi-qQ > 0/H^} = / exp -p^(^+-^) -^e:q5{-£^/2a ll (^)dp (4.53) 

0 I 2a, 2a ) o ' ’a 



Letting 



2a: 



2a( 



2a' 



2 2 
a + a. 
o 1 

■; 2 - 2 



(4.54) 



SO that 



z 2 2 ,/ z z \ 



(4.55) 



we have 



2 2 2 
“ 0 0 /2a 

Pfqj">5(j > O/H^) = / exp{-p /2o.p) -J ® °I 

“ °o °T 



Orr, a / 



2 2 2 2 2 
a^ -£ /2a^ ay2a^ 






2 2 2 
-{p^+a:)/2at 



pa 



^o< 



)dp 



(4.56) 



where 



2 , 2 

a = ea^/a 
o T o 



(4.57) 



Now the integral itself yields 1, since it is the integral 
of a probability density function. Therefore 
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P{q^-q^ >0/H^} = 



^exp 



•5 2 2.2 
V<^o 
2 2 



^T 

^exp 

<^0 



( 


r 21 


' 


1 2 


Pm 




£ 


1 - ^ 




] o 2 


^ 2 




/ 2a^ 


a 




I _ o 


o_ 





(4.58) 



Similarly for Eq. 4.51 



P{q3_^Q<0/Hj_} = / 



p „ -„ aV )/ 20 ? / q ,® 



r 



0 \ rr^ 



u(q^)dq^ 



P^(p/H^)dp 



(4.59) 



Observing first, the quantity in brackets can be expressed as 



/ 



2 2 2 AT 

P -((AT/2)^-Kj^)/2a^ /qi(Y> 



Iq — 2— u(qi)dqi. 



On 



= 1 - / 



,,AT,2 1,^ 2 , ,AT. 

% -((t) ^/2a^ ^ /qi^T' 



~i: ® 

p 






(4.60) 



Letting 



a 



^ = AT/2a^ 



(4.61) 



and making a change of variable 



X = ^i/^i 



(4.62) 
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we ha_ve that Eq. 4.60 becomes 



1 - 



p/a. 



/ 2 2,^ 

“ -{x +a,)2 

/ X e (a^x) u (x) dx 



1 - Qia^p/Oj) 



(4 



where Q(’,*) is the well-known Marcum Q function [Ref. 10]. 
Therefore Eq. 4.59 becomes 



p^qi-q^ = / 






2 /-, 2 

0 



u(p)dp 



0 '^l <^1 



dp 



(4 



From Reference 11, the integral of Eq. 4.64 becomes 






0^ ri '“ll 



dp = 



2 2 
'^o^'^l 



1 - Q 0, 




(AT/2) 

2 2 
'^o ^ '^1 



2 2 
'^o^'^l 



(V: 



. -li, 0 



2 2 ' 
'^o ^ ^1 



2 2 
a + a, L 
o 1 



T ) (AT/2) 

1 - exp - y 






+ 



2 2 



= 1 - 



2 

<^0+^1 



2 exp 



(AT/2) 



(4 



.63) 



. 64) 



.65) 
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Thus Eq. 4.59 becomes 



P{qi-q^ <0/Hi> = 2 ^ 2 ^ 



2 2 

G + 0 ’-, 
O 1 



(AT/2) 

2 (o^«J) 



(4.66) 



Now using Eq. 4.58 and Eq. 4.59 in Eq. 4.39, we have that 



P 



1 ^T 

2 — exp 



(AT/2) 



2a‘ 



1 - 



2 "| 

2 



+ 1-^ 
2 2 2 
o +o, 
o 1 



exp <- 



(AT/2) 



Recalling that £ = AT/2, using Eq. 4.55, we have 



2 



Hal +o^) 



(4.67) 



From Eq. 4.67 it is clear that in order to minimize , we 

2 2 

must maximize a +a, subject to some constraint on the colored 

o 1 

noise power. By Eq. 4.20 



2 2 2 2 2 2 

a+a, = a+ a +a+a, 
o 1 w c,o w c,l 



^22 2 
2a + a + a , 
w c , o c , 1 



(4.68 



where 



62 



2 2 

a +a T = // K (t-T)Cos oo t+Cos u t dt dx 

c,o c,l 00 ^ o o 



T T 

+ // K (t-x)Cos 00 , t Cos oo, T dt dx 
00 ^ ^ 



(4.69) 



As an example, consider the case where the power spectral 
density of the jammer is 



S (oj) — ttP [ 6 ((o~oo.) +5 (oo+oj • ) ] 
c CO] 



Under this assumption it has been shown in Appendix C, that 
Eq. 4.69 becomes 

2 2 
a +c , 
c,o C,1 



P T 
c 



L\ 



- 0 Oj)T /2 



-0Oj)T/2 



+ 



Sin(o3^-(Oj)T/2 \ 2I 

(tOj^-ojJT/2 j J 



(4.70) 



where (Oj and P^ are the 
waveform, respectively. 
Appendix C that Eq. 4.70 



frequency and the power of the jamming 
It has also been demonstrated in 
is maximum at 00 . =00 ’ or u . = to, , so 

DO ] 1 



that 




max 



P 

c 

4 



'Sin ((0^-03^) T/2 
, ( 00 ^- 00 ^) T/2 




(4.71) 



and Eq . 4.68 now becomes 
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2 



N T 
o 



2 



(4.72) 



2 2 

a + a, = 
o 1 



P 

c 



-\ 



Sin(oo^-uj^)T/2 
(03^-00^) T/2 



1 



Thus Eq. 4.67 becomes 



P 



1 

2 expi 



2 2 

-A T /4 



N T + 
o 



P T 
c 



.2 



[ 



1 + 



Sin(w^-w^)T/ 2 \ 2 




2 

A T/4N^ 

, , "/r, . I Sin(V“o>^/2Vl 

2N I -^ \ (00, -CO )T/2 / 

O «- ' i o -* 



n _ 1 - SNR 

^ Sin(oo^-co^)T /2 ^^2 

(oo^-oo^)t/ 2 j 

2 2 

where SNR = (A T/2)/N and JSR = P T/(A T/2), represent signal 

o c 

to noise ratio and jamming to signal ratio respectively. 

Observe that with JSR = 0, Eq.4.73 becomes identical to 
Eq. 4.10. This result is appealing because for the case of 
no jamming, the receiver performance should be identical to that 
of a receiver operating in white Gaussian noise interference 
only . 






2 + JSR- SNR 1 + 
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C. RECEIVER PERFORMANCE IN THE PRESENCE OF WHITE GAUSSIAN 

NOISE UNDER SINGLE CHANNEL OPERATION 

In Section B, we have analyzed the performance of the 
quadrature receiver in the presence of white and colored 
Gaussian noise. Results were specifically obtained when the 
colored noise interference was a single frequency jammer. Sup- 
pose now that the quadrature r-eceiver experiences a single 
frequency interference which corresponds to one of the signal 
frequencies, say o)^. Since the receiver makes binary decisions 
based on whether q^^ > q^ or vice versa, the presence of the 
interference at frequency will cause q^ to be greater than 
qj^ most of the time creating decision errors nearly 50% of the 
time . 

In order to prevent this type of situation from arising , 
the receiver can turn off the affected channel, or equivalently, 
make decisions based only on the output of the other channel, 
that is, based only on the size of q^^. In this section the 
performance of the quadrature receiver is analyzed assuming 
white Gaussian noise only interference, and that decisions based 
on only one channel output are made. 

Assuming that the receiver bases decisions only on the 
size of q^, the decision rule now becomes 

»1 

qj, < T (4.74) 

H 

o 

Recall from Eq . 4.6 that. 
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(4.6) 



r T -|2 r T 12 

= / r(t)Sin coj^tdt J ’*’[_/ r(t)Cos uij^tdt J 



k = 0,1 



A ,.2 . 2 



= x^ + y: 



(4 



The probability of error is 



?e = P{qi >Y/H^}P{H^} + P{q^ <Y/Hj_}P{H^} 



(4 



and assuming that P{H^} = P{Hj^} = 1/2 then Eq . 4.76 becomes 






(4 



The information bearing signals are 



/Ej_ S^(t) = A sin (cOj^t + 4)^ ) i = 0,1 

0 < t < T 



(4 



and in Section B we found that 



P(qi/Ho; 



2 ,o 2 

qq -q2_/2oj_ 



U (q. 



(4 



and 



P(qj_/H^; 



2 2 2 

q, - (£ +q, ) /2a q., e 

4e 1 u(q^: 



a. 



(4 



.75) 



.76) 



. 77) 



.4) 



.46) 



.47) 
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where e = AT/2 and 



o 



2 

1 



2 



w 



N T 
o 



4 



Thus, from Eqs . 4.19, 4.46, and 4.47, we obtain 



P 



e 



1 

2 ^ ® 

Y 



u(q^) dq^ 




' T 



2 2 2 
-(£ ^^)/2of 

e 




u(q^)dq^ 



Observe however that a threshold of y must now be defined. 

Clearly, a threshold that minimizes P should be chosen. 

e 

can be done by solving dP^/dy = 0. 

Thus 



dP 
6 

dy 



[■ 



2,„ 2 
-y /2o 

e u(y) + 



-(e^+y^)/2af 



^)U(V)] 



2^2 

2 -i e u(y) 

°1 



2 2 

-l.e 

Oi^ 



so that solution of dP^/dy = 0 yields an implicit solution 
y, namely, 



(4.21) 



(4.78) 



This 



(4.79) 



for 
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e 



(4.80) 



(t) 



2 2 
e /2o^ 



Suppose now that is the solution of Eq. 4.80 for a given 

2 

value of £ and Then 



p = I 

^o ^1 



u(q^)dq^ 



1 

^2 



/ 



Y 2 2 2 

o q^_ -(£ +q^)/2a^ /'^1^\ 



^1 



o 



u(q^) dq^ 



(4 



Letting y = q^/o^, Eq. 4.81 becomes 



P 



1 

2 ^ 



- ^/2a^ , 

Q> 1,1 






+ j 1 - 






^ 2. 

» -2 ( +y ) 

/ye I, 



(y(^))dy] 



11 / V^^i If,,/ / , 

2+je - 2 Q<=/Oi-V°l’ 



(4 



Observe that 



2a. 



(AT/2) • 

2N T/4 
o 



A^T 

2N 

o 



= SNR 



(4 



so that defining 



.81) 



.82) 



. 83) 
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(4.84) 






'<0^^ 



we have 






= Yth/2SNE 



(4.85) 



SO that the threshol(d setting equation (Eq. 4.80) becomes 



I^(y^j^( 2SNR) ) = e 



SNR 



(4.86) 



an(i Eq. 4.82 simplifies to 



= i+iexp{- I Yto^^SNR) } - jQ (/2Sfm,Y^(/2SNR)) (4.87) 

The receiver performance indicated by Eq. 4.87 is compared to 
that of an incoherent BFSK receiver that utilizes both channels 
for its decisions. (See. Eq. 4.10.) The result of this 
comparison is presented in Chapter V. 

D. RECEIVER PERFORMANCE IN THE PRESENCE OF COLORED NOISE 
UNDER SINGLE CHANNEL OPERATION 

In this section, we analyze the performance of the quadra- 
ture receiver under the assumption of single channel operation, 
as described in the previous section. Here however, it is 
additionally assumed that a jamming signal is present, whose 
energy is concentrated around the frequency (Observe that 



69 



the channel whose output is has a passband around Thus 

a jammer concentrating its energy around would significantly 

affect the output q^. Consequently, turning off or ignoring 

q would make sense under these circumstances. Hence, the 
o 

single channel operation being considered here. ) 

Our decision rule continues to be 



> 



(4.74) 




and 



P 



e 



|p { q^ > y/H^} + ip { q^ < y/H^} 



(4.77) 



Observe that due to the presence of a jammer 





(4.28) 



where 



2 

a 



T T 

// K (t-T)Sin (jo-, t Sin (jo, xdt dx 
00 



(4.20) 



As shown in Appendix C, 



F T 
c 



2 




2 



(C.12) 



4 
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when the jammer is concentrated at frequency Wj . With 

oj . = oj , Eq. C.12 becomes 
3 o' ^ 



2 



a 



1 




Sin (to^-oaQ^) T/2 
(03^-0)^) T/2 




(4.88) 



so that the probability of error is 



P 

e 




I 



1 ,,2 
2 'IH 





(4.89) 



Observe that 



o,,+o , 

w c,l 



A^T/2N 



P T 

11 c 

2 2 “2 

A T/2 



A^T/2 

N 



Sin(w -OJ 
o 

(03 -a; J T/2 

O 1 



^l)T/2l 

t72~ J 



2SNR 



1 + JSR-SI4R 






Sin(o3 -o 3,)T/21 2 
o 1 

( 03^-03^) T/2 



Defining 



SSQ = 



Sin ( ) T/ 2 

(03^-03^ )T/2 
o i 



(4.90) 



(4.91) 



we have 



p 



e 




1+JSR-SNR-SSQ 



2SNR 



1 

2 




1+JSR-SNR-SSQ 



2SNR 



4 



1-hJSR-SNR-SSQ 



2SNR 



(4.92) 



Observe that with JSR = 0, Eq.__4.92 becomes identical to Eq . 
4.87, as must be the case. 

Furthermore if the frequency separation (^^“coj^) is such 
that T/2 >> 1 or (to^-co^)T/2 = mr , where m is an integer 

then, SSQ becomes very small or zero so that the effect 
of the presence of the jamming is negligible. The numerical 
results obtained from Eq. 4.92 are very similar to those 
obtained from Eq. 4.87 as demonstrated in greater detail in 
Chapter V. 

Recall that the threshold is obtained from the solution of 
Eq. 4.86, namely 



However if our goal is to set a threshold that minimizes , 
for the case being considered here, we can solve for an opti- 
mum threshold setting by minimizing Eq. 4.92 with respect to 



Io( ) 



SNR 



(4.86) 



e 



If this procedure is carried out, we obtain the threshold 
TH 



setting equation 




(4.93) 



72 



While this result is intuitively appealing, a practical 
problem arises in that in most cases, the receiver does not 
know the operating JSR value, hence a threshold could not be 
set. 

Fortunately, computer evaluations carried out using both 
Eq. 4.86 and Eq. 4.93 to set the threshold have demonstrated 
that the resulting with thresholds set by Eqs. 4.86 and 
4.93 are almost (and for all practical purposes) identical. 



73 



V. GRAPHICAL RESULTS 



A. GRAPHICAL RESULTS FOR COLORED NOISE INTERFERENCE IN 

COHERENT M-ARY FREQUENCY SHIFT KEYED MODULATED SYSTEMS 

In Chapter III, the performance of the MFSK receiver in 
the presence of white and colored noise was derived. This 
mathematical result is used now to evaluate and graphically 
display receiver performance under the presence of white 
noise only and under the presence of- white and colored noise 
interference . 

Results are presented sequentially for values of M = 2, 4, 

8, and 16 on the performance of the M-ary FSK receiver for 
white noise as the only source of interference as well as for 
various conditions of colored noise powers 'in addition to the 
normally present WGN interference. The performance results 
for the M-ary FSK receiver presented in this section in terms 
of the probability of error are shown as the SNR changes. for 
specified values of JSR. Some representative results are 
summarized in Tables 5.1, 5.2, 5.3 and 5.4. Figures 5.1 through 
5.4 include the performance of the M-ary FSK receiver when the 
transmitted signal is interfered by white noise only, namely, 

JSR = 0. This makes it possible to evaluate the effect of the 
jamming on the receiver in comparison to the case in which 
WGN is the only source of interference. These results have 
been obtained by evaluating Eq. 3.47. 
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TABLE 5.1 



PERFORMANCE OF 2-FSK RECEIVER 



The Receiver 


P 

e 

SNR (DB) 


-10.0 


-5.0 


O 

• 

o 


5.0 


10.0 


JSR = 0 


0.3759 


0.28 6-9 


0.1586 


0 .0376 


0 . 0008 


JSR = 0 db 


0.3815 


0 .3120 


0.2397 


0.1917 


0 . 1702 


JSR = 5 db 


0 . 3914 


0 . 3454 


0.3120 


0.2959 


0.2899 


JSR = 10 db 


0.4115 


0 . 3914 


0.3815 


0.3778 


0 . 3765 


JSR = 15 db 


0.4384 


0 .4327 


0.4305 


0.4298 


0 .4295 
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TABLE 5 . 2 



PERFORMANCE OF 4-FSK RECEIVER 



The Receiver 


P 

e 

SNR (DB) 


-10.0 


-5.0 


0.0 


5.0 


10.0 


JSR = 0 


0.6223 


0.5132 


0.3222 


0.0915 


0.0022 


JSR = 0 db 


0.6262 


0.5313 


0.'3995 


0.2804 


0.2153 


JSR = 5 db 


0.6326 


0.5598 


0.4861 


0.4395 


0.4194 


JSR = 10 db 


0.6478 


0 . 6082 


0.5825 


0.5712 


0.5671 


JSR = 15 db 


0.6734 


0.6598 


0.6538 


0 . 6517 


0.6510 
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TABLE 5 . 3 



PERFORMANCE OF 8-FSK RECEIVER 



THE RECEIVER 




SNR 


P 

?DB) 






o 

o 
1 — 1 
1 


-5.0 


0.0 


5.0 


10.0 


JSR = 0 


0.7778 


0.6794 


0.4755 


0.1617 


0.0048 


JSR = 0 db 


0.7792 


0.6885 


0.5261 


0.3246 


0.1885 


JSR = 5 db 


0.7820 


0.7047 


0.5958 


0.4992 


0.4471 


JSR = 10 db 


0.7894 


0.7384 


0.6914 


0.6648 


0.6540 


JSR = 15 db 


0.8056 


0.7834 


0 . 7709 


0.7658 


0.7641 
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TABLE 5.4 



PERFORMANCE OF 16-FSK RECEIVER 



THE RECEIVER 


P 

SNR (DB) 


-10 . 0 


-5.0 


o 

o 


o 

• 

IT) 


10.0 


JSR = 0 


0.8715 


0.7949 


0.6083 


0.2455 


0 . 0093 


JSR = 0 db 


0.8720 


0.7987 


0.6354 


0.3621 


0.1374 


JSR = 5 db 


0.8731 


0.8062 


0.6796 


0.5174 


0 . 4019 


JSR = 10 db 


0.8763 


0.8243 


0.7556 


0.7000 


0 . 6733 


JSR = 15 db 


0 . 8839 


0.8553 


0.8329 


0.8217 


0.8175 
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MFSK (M=2) 




Figure 5.1. Performance of M-ary FSK for M = 2 
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MFSK (M=4) 




-10.0 0.0 10.0 20.0 ao.o 40.0 



SNR DB 



Figure 5.2. Performance of M-ary FSK for M = 4 
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MFSK (M=8) 




Figure 5.3. Performance of H-ary FSK for h = 8 
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MFSK (M= 16 ) 




I I I 

• 10.0 0.0 10.0 20.0 30.0 40.0 



SNR DB 



Figure 3.4. Performance of M-ary FSK for M = 16 
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B. GRAPHICAL RESULTS FOR NON-COHERENT BINARY FREQUENCY 
SHIFT KEYED SIGNAL DETECTION IN THE PRESENCE OF 
COLORED NOISE 

In Chapter IV, the performance of the quadrature receiver 
operating in the presence of white and colored noise was 
derived. The mathematical results are now used to evaluate 
and graphically display receiver performance under various 
conditions of signal and noise powers. 

First, results are presented for the case in which white 
noise is the only source of interference. This yields the 
well-known probability of error curves for the standard quadra- 
ture receiver for non-coherent BFSK . These are presented in 
Fig. 5.5, along with a corresponding plot of the probability of 
error of the quadrtature receiver in which only one channel 
output is used to make binary decisions. 

Additionally, the performance of the quadrature receiver 
operating in the presence of white and colored noise is evalu- 
ated under dual channel and single channel operation. Under 
single channel operation, it is assumed that the colored noise 
jamming concentrates its energy around one of the FSK operating 
frequencies, and that the receiver is able to make a determinis- 
tic as to which "channel is being jammed" so that the outputs 
of this channel are ignored in the process of making decisions. 
Evaluations are carried out using receiver thresholds that are 
dependent as well as independent of jamming power levels. 

(Both cases are considered separately.) The performance of 
the quadrature receiver in the presence of noise and the jamming 
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QUADRATURE RECEIVER JSR=0 




Figure j.d. Performance of the quadracure receiver for JSR = 0 
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waveform described in this section in terms of the probability 
of error is calculated as the SNR changes for specified values 
of JSR. Some important results are summarized in Table 5.5 
for JSR = 0 and in Tables 5.6-5.10 as JSR takes on values 
of 0.0 db, 5.0 db , 10.0 db, 15.0 and 20.0 db, respectively. 

In Figure 5.5 the per formance of the standard quadrature receiver 
and the single channel operation of the quadrature receiver is 
plotted when the transmitted signal is interfered by white 
noise only. The theoretical performance of the standard 
quadrature receiver is calculated from Equation 4.10, and the 
performance of the quadrature receiver under single channel 
operation is calculated from Equation 4.87. 

In Figures 5.6-5.10, the performance of the standard 
quadrature receiver and the quadrature receiver under single 
channel operation with the threshold dependent as v/ell as 
independent of the jamming power level is plotted when the 
transmitted signal is interfered by white noise and by the 
jamming waveform having Power Spectral Desnity given by Equation 
C.7. Each of the figures corresponds to a specific value of 
JSR as shown in the headings. The performance of the standard 
quadrature receiver is calculated from Equation 4.73. The 
theoretical results for the single channel operation of the 
quadrature receiver with a threshold that is independent of the 
jamming power (Eq. 4.86) is calculated from Equation 4.82, and 
Equation 4.92 is used to compute performance of the same re- 
ceiver when the threshold is dependent on the jamming power 
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TABLE 5.5 



PERFORMANCE OF THE QUADRATURE RECEIVER JSR = 0 



THE RECEIVER 


P 

e 

SNR (DB) 




-10.0 


-5 . 0 


0.0 


5.0 


10.0 


15.0 


Standard 

Operation 


0.4756 


0.4268 -■ 


0. 3032 


0. 1028 


0.0033 


0 . 0000000b 


Single 

Channel 

Operation 


0 .4820 


0.4460 


0.3531 


0 . 1806 


0.0268 


0.00009 



TABLE 5 . 6 

PERFORMANCE OF THE QUADRATURE RECEIVER JSR = 0 DB 



THE RECEIVER 


P 

e 

SNR (DB) 




1 

O 

o 


-5 . 0 


0 . 0 


5 . 0 


10 . 0 


15.0 


Standard 

Operation 


0.4767 


0.4361 


0. 3582 


0.2709 


0. 2172 


0. 1952 


Single 

Channel 

Operation 


0 .4820 


0.4460 


0.3531 


0 . 1806 


0 . 0268 


0.00009 
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TABLE 5 . 7 



PERFORMANCE OF THE QUADRATURE RECEIVER JSR = 5 DB 



THE RECEIVER 






P 

e 

SNR (DB) 








-10.0 


-5.0 


o 

o 


5.0 


10 . 0 


15.0 


Standard 

Operation 


0.4788 


0.4499 


0.4119 


0.3841 


0.3713 


0.3667 


Single 

Channel 

Operation 


0.4820 


0.4460 


0.3531 


0.1806 


0.0268 


0.00009 



TABLE 5.8 

PERFORMANCE OF THE QUADRATURE RECEIVER JSR = 10 DB 



THE RECEIVER 


P 

e 

SNR (DB) 




1 

O 

O 


-5.0 


0.0 


5.0 


10 . 0 


15.0 


Standard 

Operation 


0.4836 


0 . 4702 


0.4600 


0.4551 


0.4533 


0.4527 


Single 

Channel 

Operation 


0.4820 


0.4460 


0. 3531 


0.1806 


0.0268 


0.00009 
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TABLE 5.9 



PERFORMANCE OF THE QUADRATURE RECEIVER JSR = 15 DB 



THE RECEIVER 


P 

e 

SNR (DB) 




o 

o 

1 


-5.0 


0.0 


5.0 


10.0 


15.0 


Standard 

Operation 


0.4904 


0.4869 


0.4853 


0.4847 


0.4845 


0 .4844 


Single 

Channel 

Operation 


0.4820 


0.4460 


0.3531 


0.1806 


0.0268 


0.00009 



TABLE 5.10 

PERFORMANCE OF THE QUADRATURE RECEIVER JSR = 20 DB 



THE RECEIVER 


P 

e 

SNR (DB) 




-10.0 


-5.0 


0.0 


5.0 


10.0 


15.0 


Standard 

Operation 


0.4958 


0.4953 


0.4951 


0.4950 


0.4950 


0 . 4950 


Single 

Channel 

Operation 


0.4820 


0.4460 


0.3531 


0. 1806 


0 . 0268 


0.00009 
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QUADRATURE RECEIVER JSR=0 DB 




- 10.0 0.0 10.0 20.0 30.0 40.0 



SNR DB 



Figure 5.6. Performance of the quadrature receiver for 
JSR = 0 db 
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Figure 5 



QUADRATURE RECEIVER JSR=5 DB 




.7. Performance of the quadrature receiver for JSR = 



5 db 



90 



QUADRATURE RECEIVER JSR=10 DB 




SNR DB 



Figure 5.8. Performance of the quadrature receiver for JSR = 10 do 
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QUADRATURE RECEIVER JSR=15 DB 




SNR DB 



igure 5.9. Performance of Lhe quadrature receiver for JSR = 



15 db 



QUADRATURE RECEIVER JSR=20 DB 




Figure 3.10. Performance of the quadrature receiver for JSR = 20 db 
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level. As pointed out in Section D of Chapter IV, the proba- 
bility of error calculated from Equation 4.92 with the thres- 
hold set by Equations 4.86 and 4.93 show almost identical 
results . 

Tables 5.5 through 5.10 demonstrate that the performance 

of the quadrature receiver under single channel operation is 

unaffected by changing values of JSR. This is due to the fact 

that for CO / co, and T values used in the simulation, the value 
o 1 

of SSQ term in Eq. 4.92 is identical to zero. Thus in order 
to demonstrate the effect of the jammer on the receiver under 
single channel operation, the value of the jamming frequency 
coj has been allowed to vary from to^ all the way up to toj^. 

Thus, in place of the SSQ term as defined in Eq. 4.92, we use 
the modified term 

Sin (co . - 0 ), ) T/2 

SSQ = — ^ — rrmTo ^ < co . < co, 

(c0j-co^)T/2 o - 3 — 1 

The results of these modifications are presented in Fig. 5.11 

and Fig. 5.12 where the probability of error of the receiver 

given by Eq. 4.93 is evaluated for JSR = 5 db and JSR = 10 db, 

respectively, where the jamming frequency (co^) is allowed to 

take on values co . = co (which corresponds to the results given 

30 

by Eqs. 4.91 and 4.92 without modification), and values of 
Uj = 3 (coj^ + 00 ^) /4 and cOj = co^^. Some of the important results 
obtained are summarized in Tables 5.11 and 5.12 for JSR = 5 db 
and JSR = 10 db respectively. 
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B F S K (JSR=5 DB) 




*^- 10.0 0.0 10.0 20.0 30.0 40.0 



SNR DB 



Figure 5.11. Perforraance of the quadrature receiver sing 
operation for dizferent jaroining frequencies 
JSR = 5 db 



e channel 
and 
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B F S K (JSR=10 DB) 



Figure 5.12 




Performance of the quadrature receiver single channel 
operation for different jamming frequencies and 
JSR = 10 db 

OG 




TABLE 5.11 



PERFORMANCE OF THE QUADRATURE RECEIVER SINGLE CHANNEL OPERATION 
FOR DIFFERENT JAMMING FREQUENCIES AND JSR = 5 DB 





SNR DB 


THE RECEIVER 


-10.0 


0.0 


10.0 


20.0 


30.0 


0) . = 0) 

D o 


0.4820 


0.3530 


0.0268 


1 X lo"^ 


-9 

1 xio 


0). = 


0.4821 


0.3589 


0.0676 


0.0049 


0.0023 


U3j = 0)^ 


0.4853 


0.4476 


0.4297 


0.4272 


0.4269 



TABLE 5.12 

PERFORMANCE OF THE QUADRATURE RECEIVER SINCLE CHANNEL OPERATIC.’. 
FOR DIFFERENT JAMMING FREQUENCIES AND JSR = 10 DB 





SNR DB 


TliE RECEIVER 


-10.0 


0.0 


10.0 


20.0 


30.0 


e 

II 


0.4820 


0.3531 


0.0268 


-9 

1 X 10 


1 X lo”^ 


rjo . = ) 

j 4 1 o 


0.4823 


0.3699 


0.1491 


0.0849 


0.0773 


03 . = COt 

J 1 


0.4894 


0.4781 


0.4755 


0.4752 


0.4751 
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VI. CONCLUSIONS 



The analysis carried out in this thesis presents the 
application of concepts derived in statistical communication 
theory, specifically in the theory of signal detection under 
the assumption of colored noise interference. The performance 
of digital receivers in terms of probability of error is 
determined when the receivers operate in the presence of white 
and colored Gaussian noise. Three techniques are examined 
separately, one for MPSK modulation, another for coherent 
MFSK modulation and the last one for (incoherent) BFSK 
modulation. 

The mathematical model of the jamming waveform proposed, 
consists of colored Gaussian noise of different spectral shapes 
and power content. 

For MPSK modulation, a mathematical result on the performance 
of the (coherent) receiver in the presence of WGN and colored 
noise jamming was derived. The complexity of the result along 
with the many possible trade-offs involving spectral shapes, 
power levels and frequencies of operation made it impossible 
to address in this thesis the issue of optimum jamming strate- 
gies for MPSK. 

For MFSK modulation results on the effect of the coherent 
receiver, were derived. A simple assumption was m.ade on the 
spectrum of- the jamming. By assuming that each signal frequency 
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was interfered with a tone subject to a total jamming power 

constraint, the receiver was evaluated for different values 

of SNR, JSR, and M. The results demonstrate that this 

form of jamming can be quite effective or that significant 

increases on can be achieved even at low JSR values, 
e 

For the case of BFSK modulation, the quadrature receiver 
was analyzed under two conditions of operation, standard 
operation and single channel operation, in the presence of 
colored noise jamming with different power levels. The single 
channel operation was introduced as a method for mitigating the 
effect of a single tone jammer at one of the carrier frequen- 
cies. When no jamming is present, single channel operation 
performs slightly worse than standard receiver (both channels) 
operation. However, in the presence of jamming, single channel 
operation is superior to standard operation because the receiver 
is capable of eliminating much of the jammer energy and its 
effect by -ignoring the output of the jammed channel during 
single channel operation. As pointed out in Chapter IV, the 
effect of the jamming waveform on the receiver under single 
channel operation depends strongly on the jamming frequency 
chosen. For the single channel operation, it was assumed that 
the jamming is present at one of the two signal frequencies, and 
that the receiver turns off the channel affected. Thus, deci- 
sions are made based only on the output of the unaffected 
channel. However, if under this condition of operation the 
jamming changes its frequency ojj in such a way as to "move 
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closer" to the frequency of the unaffected channel, it has 
been demonstrated that the receiver probability of error in- 
creases as o)j approaches the frequency of the unaffected 
channel . 
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APPENDIX A 



DETAILED INVESTIGATION OF THE VARIANCES OF Vp_ AND V; 
CONDITIONED ON HYPOTHESES H j 



Let 






and 



(t) 



0<t<T £=1,2 



otherwise 



T 



F{t|;^(t) } 1=1,2 



Thus 



OO 00 



c,£ 



// K (t-T)i|;'J(t)ii;J(T)dt dx 



— oo-OC 



oo oo 



/ / ^S "'■^dwi|j'^(t)ii^”^(T)dt dx 



— OO — oo 



'H' ' 



^ / S^ (x)i|;'J(-u})i|;'J(aj)daj 



where K (x) 
c 



Sc(w 



Now 



(A. 1) 



(A. 2) 



(A. 3) 
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i[/2f 



-2irr/T)T/2 ^ ) T/2 

(0)-^)T/2 



+ /2T e‘^("^2TTiVT)T/2 S^(^ +^?')T/2 

(w+^)T/2 



(A. 4) 



and 



i|^2 (w) 



-jcoT/4n ,T, . 



(A. 5) 



T T 

Because of the relationship between and i |^2 (w) , it is 

clear that 



T T 

(-co) (go) 



T T 

ip2 ^2 



(A. 6) 



so that 



2 2 A 2 

y , = a - = a 

c , 1 c , 2 c 



(A. 7) 



Thus indeed, a = c. Observe also by similiar arguments, that 



. VJU , 

?1,2 = ^ / S ^ (o)).l;^(-w)iiJ2(oa)da) = ^ f (co)i|;^(-a))e^‘^‘'^%^(co)daj 



oo 2 

" ^ ^ Cos ^ doj 



(A. 8: 
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where 



ki(aj) 1^ 



psin(^^|-nT) -|2 [- sin(^ + nT)1 2 

2 ,toT ,ooT^ ; I 

L {— -nT) J L (— + riT) J 



Sin (^- nr) 
(^-nr) 



Sin(^ + nT) 

(^ + nV) 



(A. 9) 



So it is clear that in general. 



0i 2 will not be zero. 
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APPENDIX B 



DETAILED INVESTIGATION OF THE BEHAVIOR OF THE PRODUCTS OF 

(-0)) AND S^(u) 



We have defined 



F{S^ (t)p(t) } 



r y / \ 

-t 5 “LL ((jo-o) -go . ) 

Z ^ J 



L {(jO +03 +(jO . ) ] 
^ J 



(B.l) 



Then 



F{S! (t)p(t) } X F{Sj;^(t)p(t) } 



= (-co)Sj;. (co) 



AT AT 

— [L(-o)-a) “CO ■ ) + L ( -co+co^+co . ) ] [L { (jo-co “co, ) + L { cu+co^+cot_ ) 

Z CJ CjZ CK C K 



AT 

(— ::p) [L(- 50 -C 0 ~(jO . ) L ( (jO”CO *C0, ) + L ( — CO + CO +CO.)L(cO — CO “COi ) 

2. Cj CK CJ CK 



+ L ( ~ CO ~co ~u) . ) L (co+cl) +00 . ) + L (~oo+co +co • ) L (co+co +co, ) 

CJ CJ CJ CK 



(B.2) 



Observe that for reasonably large values of oo^ , the first and 
the last term in this expression vanish, and we are left with 
the products 
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j-(o)|^-(jjj) T/2 Sin(o)-u)^-(j0j)T/2 Sin(ur-aj^-Wj^)T/2 

( CO-CO -u}.)T/2 
^ 3 



((o-co^-a^)T/2 



+ e 



■j (c0j^-C0j)T/2 Sin((of(o^+(Oj)T/2 Sin(cof(o^+a:^)T/2 

((of(o^+(Oj)T/2 ((ofco^+(Oj^)T/2 



(B. 3) 



Now focusing on the first teritr of Eq. B.3, which has significant 

components for co in the neighborhood of co^ , we see that if 

j-k >> 1 then there is essentially no overlap between sine 

functions. Therefore the product S 1 (-co) S, (co) is zero for 

3 ^ 

j 7 ^ k- 

For k = j ±1, we have 



Sin ( co“co_~co'- ) T/2 Sin (co-co^~coi, ) T/2 

C J C K 

= ^ Cos (coj^-cOj ) T/2 - ^os (2(o-2co^-(Oj 

and when co is in the neighborhood of co^ / 
approximately 



-t0j^)T/2 (B.4) 

the product becomes 



j (‘Ok ^ r Cos (coj^-co • ) T/2 - Cos (co .+co^) T/2‘ 



CO j T/2- cOj^T/2 



j (k- j ) AcoT/2 r Cos (k- j ) AcoT/2 - Cos (k- j ) AcoT/2 



= e 



(k-^)T/2(j -^^)T/2 J 



(B.5) 
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The orthogonality condition on the signals required that 
Ao)T = IT or Aa)T/2 = 'it/2, so that Eq. B.5 becomes (approximately) 



(T/2)^ 



Eq . B . 6 
that for 
j ^ k. 



E 



Gos(±tt/2) -• Cos(2j ±1)tt/2 






(D — —) 



±j fCos jTT Cos 7t/2 isin jir Sin Tr/2~\ . 

= (T/2)H J 

is zero for all values of the integer j , so we have 
AojT = IT, the product (- 00 ) Sj^ (oj) is equal to zero for 
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APPENDIX C 



DETAILED INVESTIGATION OF THE VARIANCES 

AND , DUE TO COLORED NOISE 
c , o c , 1 



Let us define 



Pci(t) = Cos oj^t i = 0,1 0 £ t £ T (C.l) 



Then Eq. 4.69 becomes 



2^2 
a + cr , 
0,0 0,1 



fj (t-T)P (t)P (T)dt dT 



00 oo 



+ // K^(t-T)P^^(t)P^j_(x)dt dT 



— s - O - X 



3 C ^ 

j ^ jj e^^^^ "'■^P _ (t)P (x)dt dx dm 

2 tt ^ c ^ CO CO 

—CO — oo-O^ 



oo o 

= ^ / S^(u.) // 



00 oo ^ . 

+ i / S (m) // (t)P„, (x)dt dx 



2tt ^ c 

— OO — OO-'^ 



cl cl 



S I Scto)t|P^to)r + 



(C.2) 



where 
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Cos co^te 



-jtjjt 



dt 



^ j (to^-co) T/2 Sin ((jj^-co) T/2 
2 ® (co^-Gj) T/2 



T -j (co^+to) T/2 Sin (oj^+co) T/2 
^ 2 ® (co^+uj)T/2 



(C.3) 



and 



S^(oj) ^ K^(t) 

c c 



Furthermore 






= 



Sin T/2 

(o)^-a))T/2 



^ 2 ’ 



Sin (oj^+co) T/2 
( 0)^+0)) T/2 



T, 2 



(f) 



2 Cos OJ^T 



Sin (oo^-oj) T/2 
(cj^-o))T/2 



Sin (ojj^+oj) T/2 
(ojj^+oj) T/2 



(C.4) 



i = 0,1 



The third term in Eq . C.4 can be assumed for all practical 
purposes to be zero. In essence, we require that co^ >> 2 tt/T, 
i = 0,1 for the approximation to be correct. 

Consider now the case where 



S (oj) = K[6(oj-oj-) + 5 (gj+(d-)] 

c DO 



(C. 5) 



Since 
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oo 

^ / S (a))doj 

c 2 tt c 



K 

2tt 



[ 6 (oa-ooj ) 



6 (to+co . ) ] du = — 

J 'fT 



(C. 6) 



then K = ttP^ so that we use 



S (oj) = TT p [6(w-o).) + 6(aj+to.)] 

c c D D 



(C. 7) 



From Eq . C.2, we now have 



a + a 



c,l 









(C. 8) 



A ssuming that ojj will always be in the vicinity of gOq and 



co^, we can state that 



, ,,2 - 

l^r-nKH *2* ' (u -u )T/2 ' 



CO J ' 



(c. 9 ) 



o ] 






~ ,T,2 fSin(a.^-o..)T/2^j 2 

hJ \ (w^-c0j)T/2 / 



(C. 10) 



^ y .Sin(co -u).)T/2\ 2 



1^01*'“: 



(C. 11) 
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Thus 



c,o 



+ a 



c,l 






Sin (o) -to . ) T/2 \ 2 

o 1 ] 

(to -to . ) T/2 / 

o 3 ' 



T 2 

2 (|)' 



Sin (to^-tjj ) T/2 

(to^-tOj)T/2 



Sin(co^-„.)T/2 ^2 
4 (io^-aiyT/2^ / 



Sin (to^-tOj ) T/2 
(to^-to j ) T/2 




(C. 12) 



In order to maximize the quantity in brackets as a function of 
to j , we need to take derivatives of the expression and set it 
equal to zero. The result of this operation leads to a 



maximum at values of tOj = to^ or to^ = to^. Therefore, the 
maximum value becomes 



/ 2 2 \ 


Po"' r 


/ Sin (tOj^-to^) T/2 




y'^c , o '^c , 1 / max 


4 L 


, \ (to^-to^)T/2 


/ + ij 



(C. 13) 
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